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£?*sti^yOr&oyfi«sic lM»ric«liw» is • for* of f itaio-f i)» laliricitusn 
Mere elastic oef&rswtiM of Hearing surfaces Hecfi*es significant. It is 
usually associated with ftighlj stressed sechine co«(ionents of low conform- 
ity, SiiCh aS gears and ra^l5ar»g-elei^nt Hearings. This IMrication «ecfianisn> 
IS also encourterec with soft Hearing naterials, Sw^ aS ruHHer seals and 
5 tires. The coaeorj factor in ttaese applications is that local elastic ae- 
2 , formation of tne solids provides coherent fluid films aod thus asperity in- 
teraction IS prevented. 

Historically, elastohydrodynamic lubrication may He viewed as one of 
the major developments in the field of tribology in the twentieth century, 
ii. iKit only revealed the existence of a previously unsuspected regime of 
luHrication in highly stressed and nonconformai machine elements, such as 
gears and rolling-element Hearings, but it brought order to the understand- 
ing of the co^lete spectrum of lubrication ranging from boundary to hydro- 
dynasic. 

h chronology of the research on elastohydrodynamic lubrication (£HL) of 
elliptical contacts that the author has been involved with is given in 
table i. In this table the major research efforts and some of the results 
are given for the years 15#74 to iy«0. This particular lecture concentrates 
on folly flooded conjunctions, A fully flooded conjunction is one in which 
the film ihicknass is not significantly changed when the Mount of lubricant 
is increased. 

h brief description of the relevant equations used in the elastohydro- 
dynamic lubrication analysis is given in this lecture. The most important 
practical aspect of the elastohydrodynMic lubricat-on theory is the deter- 
mination of the mininium film thickness within the contact. The maintenance 
of a fluid film of aflequete magnitude is an essential feature of the correct 
operation of lubricated machine elements. The results to be presented show 
the influence of contact geometry on minimum film thickness as expressed by 
the ellipticity parameter and the dimensionless speed, load, and aiaterials 
parameters. These results are applied to materials of high elastic modulus. 
Such as metal, and to materials of low elastic modulus, such as rubber. The 
solutions for materials of high elastic modulus are sometimes referred to as 
"hard £HL," and the solution for materials of low elastic hiodulus as "soft 
£HL," In addition to che film thickness equations that are developed, con- 
tour plots of pressure and film thickness are presented that show the essen- 
tial features of elastohydrodynamical ly lubricated conjunctions. In partic- 
ular the crescent-shaped region of minimum film thickness, with its sloe 
lobes in which the separation between the solids Is a minimum, clearly 
emerges in the numerical solutions. These theoretical solutions for film 
thickness have all the essential features of previously reported experi- 
mental observations based on optical interferooietry. Correlation between 
theory and experiments is also made. 
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PQ density at atmospheric conditions 
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Film thickness equation 
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The prot)*e« is to calcuine the pressure aistripution in the contact end «t 
the se*e tiae allow for the effects that this pressure will have on the 
properties of the fluid ano on the gecmetry of ti^ elastic solids. The 
solution win also provide the shape of the lubricant fila, particularly the 
ainiMja clearance between the sulids. m detailed description of tt^ 
elasticity aodel used is given in Oowson and Haarock (ly76), and the coa- 
plete £hl theory is given in haarock and Oowson (1^76). 


OlNENSlCNoESS GRtMjPlNG 

The variables resulting froa the isothen»«l elliptical contact theory 
developed in Haiarock and Oowson {i97oj are 

£' effective elastic modulus, h/m^ 

F noraol applied load, N 

h film thickness, a 

R* effective radius in x (motion) direction, m 

Ry effective radius in y (transverse) direction, a 

u mean surface velocity in x direction, a/s 
« pres sure-vi sees ity coefficient of fluid, 

no atmospheric viscosity, N s/m^ 

From the eight variables the following five diiriens ion less groupings were 
established: 

Dimensionless film thickness 



where 



Ellipticity parameter 


where 





0.64 



Equation (y) was obtained from Brewe and Hamrock (iy77). 
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Qiaiens ion less loaJ porasieter 



where 

F ~ normal opplieJ loaO 
OlBiensionless spcea poronw-ter 


U = 



(il) 


li^) 


Diniensionless triateriols parameter 


6 = a£* (13) 

The diniensionless film thickness can thus be written as a function of 
the other four parameters 


H = f(k.U.W.G) (14) 

The influence of the dimensionless parameters k» u» W, and G on minimum 
film thickness Hq,in is presented later for both hard and soft contacts. 

The set of dimensionless groups {H,k,U,W,G} is a useful collection of 
parameters for evaluating the results presented in this lecture. It is also 
comparable to the set of dimensionless parameters used in the initial elas- 
tohydrodynamic analysis of line contacts, and it has the merit that the 
physical significance of each term is readily apparent. However, a number 
of authors, for example, Moes (1965-66) and Theyse (1966), have noted that 
this set of dimensionless groups can be reduced by one parameter without any 
loss of generality. This approach is followed in the fourth lecture, where 
the film thicknesses to be expected in each of the four regimes of fluid- 
film lubrication are presented graphically. 


HARO EHL RESULTS 

By using the numerical procedures outlined in Hamrock and Dowson (1976) 
the influence of the ellipticity parameter and the dimensionless speed, 
load, and materials parameters on minimum film thickness has been investi- 
gated for hard EHL, fully flooded contacts (Hamrock and Dowson, 1977). The 
ellipticity parameter k was varied from 1 (a ball-on-plate configuration) 
to 8 (a configuration approaching a rectangular contact). The dimensionless 
speed U was varied over a range of nearly two orders of magnitude, and the 
dimensionless load parameter W over a range of one order of magnitude. 
Situations equivalent to using solid materiols of bronze, steel, and silicon 
nitride and lubricants of paraffinic and naphthenic oils were considered in 
an investigation of the role of the dimensionless materials parameter G. 

The 34 cases used to generate the minimum-film-thickness formula are given 
in table II. In the table Hm^n corresponds to the minimum film thick- 
ness obtained from the EHL elliptical contact theory given in Hamrock and 
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OOMSon (1976). The ninlimjov-f 1 ]»-thickness fomula obtained froffi a 
least-squares fit of the data was first given in Hamrock and Oowson (1977) 
and Is given here as 


= 3.03 


U3) 


In- table 11 is the minimum film thickness obtained from equation 
(15). The percentage difference between Is expressed 
by 




100 


(16) 


In table 11 the values of 0^ are within *6 percent. 

It IS interesting to conipare the elliptical-contact, minimum-f 1 lut- 
thlckness formula (eQ (15)) with the corresponding equation generated by 
Dowson (1968) for nominal or rectangular contacts 


min,r r 


( 17 ) 


wnere 


M 


r 



( 1 «) 



load per unit length 


(19) 


The powers of U, G, and W in equation (15) and U, G, and Wp in equa- 
tion (17) are quite similar considering the different numerical procedures 
on which they are based. It is also worth noting that the power of W in 
equation (15) is extremely close to the value of -0.074 proposed by Archard 
and Cowking (1965-66) in their experimental studies of elliptical contacts. 

In practical situations there is considerable interest in the central 
as well as the minimum film thickness in elastohydrodynamic contacts. This 
iS particularly true when traction is considered since the surfaces in rela- 
tive motion are separated by a film of almost constant thickness that is 
well represented by the central value over much of the Hertzian contact 
zone. The procedure used in obtaining the central film thickness was the 
same as that used in obtaining the minimum film thickness. The central- 
film-thickness formula for hard Ehl contacts obtained from the numerical 
results is 



2.69 - 0.61 


-0 . 73k > 

e ) 


(i^O) 


A representative contour plot of dimensionless pressure is shown in 
figure 1 for k = 1.26, U = 0.168x10-^^, w = O.lllxlO'^, and G - 4522. 

In this figure and in figure 2, the + symbol indicates the center of the 
Hertzian contact zone. The dimensionless representation of the X and Y 
coordinates causes the actual Hertzian contact ellipse to be a circle re- 
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gardless of the value of the ellipticity parameter. The Hertzian contact 
circle is shown by asterisks. On the figure is a key showing the contour 
labels and each corresponding value of diitiensionless pressure. The inlet 
region is to the left and the exit region is to the right. 

The pressure gradient at the exit end of the conjunction is ntuch larger 
than that in the inlet region. In figure 1 a pressure spike is visible at 
the exit of the contact. 

Figure 2 shows contour plots of film thickness for k = 1.25» 

U * 0.168x10“^^, U = 0.111x10'^, and G = 4522. In this figure two 
minimum-f ilm-thickness regions occur in well-defined side lobes that follow* 
and are close to, the edge of the Hertzian contact ellipse. These results 
reproduce all the essential features of previously reported experimental 
observations based on optical interferometry {Cameron and Gohar, 1966). 

The variation of pressure and film thickness in the direction of roll- 
ing quite close to the X axis near the midplane of the conjunction is 
shown in figure 3 for three va^ s of U. The values of the dimensionless 
load, materials, and ellipticity parameters were held constant at k = 6, 

U = 0.737 x 10"6, and 6 * 4522. In figure 3(a) the dashed line corresponds 
to the Hertzian pressure distribution. This figure shows that the pressure 
at any location in the inlet region rises as the speed increases, a result 
that is also consistent with the elastohydrodynamic theory for line or 
rectangular contacts. Furthermore, as the speed decreases, the height of 
the pressure spike decreases and the hydrodynamic pressures gradually 
approach the semiel 1 iptical form of the Hertzian contact stresses. Note 
that the location of the pressure spike moves downstream toward the edge of 
the Hertzian contact ellipse as the speed decreases. For nominal line or 
rectangular contacts Oowson and Higginson (1966) showed results similar to 
those in figure 3(a) . 


The typical elastohydrodynamic film shape with an essentially parallel 
section in the central region is shown in figure 3(b). There is little sign 
of a reentrant region in this case, except perhaps at the lowest speed. 

Also, there is a considerable change in the film thickness as the dimension- 
less speed is changed, as indicated by equation (15). This illustrates most 
clearly the dominant effect of the dimensionless speed parameter u on the 
minimum film thickness in elastohydrodynamic contacts. 

The variation of pressure and film thickness in the direction of motion 
along a line close to the midplane of the conjunction is shown in figure 4 
for three values of dimensionless load parameter. The values of the dimen- 
sionless speed, materials, and ellipticity parameters were held fixed at 
U = 0.168x10'^^, G = 4522, and k = 6. Note from figure 4(a) that the 
pressure at any location in the inlet region falls as the load increases. 

For the highest load (W = 1.106x10-6), figure 4(b), the film thickness 
rises between the central region and the outlet restriction. This reentrant 
effect is attributed to lubricant compressibility. Note also that at 
W « 0.5228x10-6 the film thickness is slightly smaller than at 
W » 1,06x10-6. This somewhat curious result is linked to the fact that 
the location of the minimum film thickness also changes drostically over 
this load range. At the lower load the minimum film thickness is located on 
the midplane of the conjunction downstream from the center of the contact; 
at the higher load it moves to the side lobes as described earlier. 
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COMPARISON BETWEEN THEORETICAL AND EXPERIMENTAL FILM THICKNESSES 


The minimum- and central-film-thickness formulas developed in the 
previous section for fully flooded conjunctions are not only useful for de- 
sign purposes, but they also provide a convenient means of assessing the 
influence of various parameters on the elastohydrodynamic film thickness. 

For the purpose of comparing theoretical fi'm thicknesses with those found 
in actual elastohydrodynamic contacts, table HI was obtained from Kunz and 
Winer (1977). The experimental apparatus consisted of a steel ball rolling 
and sliding on a sapphire plate; this generated a circular conjunction, or 
an ellipticity parameter of unity. Measurements were made by using the 
technique of optical interferometry. Table HI shows the results of both 
calculations and measurements for three lubricants in pure sliding, each 
under two different loads and three speeds. The notation and 

He is used to denote the dimensionless film thickness calculated from 
equations (15) and (20), respectively. The measured minimum and central 
film thicknesses obtained from Kunz and Winer (1977) are denoted by 
and Hj-. Figures 5 and 6 show comparisons between' the calculated and mea- 
sured film thicknesses for the two loads shown in table HI. 

For the smaller load (W = 0.1238x10-6) the results shown in figure 5 
compare remarkably well if we bear in mind the difficulties associated with 
the experimental determination of such small quantities under arduous 
conditions and the error assc*iated with the complex numerical evaluations 
of elastohydrodynamic conjunctions. The ratios between the central and 
minimum film thicknesses are similar for the calculations and the measure- 
ments, and the dependence of film thickness on speed thus appears to be well 
represented by equations (15) and (20). 

For the larger load (W = 0.9287x10-6) the agreement shown in figure 6 
is not quite so good, with the theoretical predictions of film thickness 
being consistently larger than the measured values. This discrepancy is 
sometimes attriouted to viscous heating, as discussed by Greenwood and 
Kauzlarich (1973), or to non-Newtonian behavior, as discussed by Moore 
(1973). Viscous heating appears to enjoy the most support. Since the mea- 
surements were made in a condition of pure sliding, thermal effects might 
well be significant, particulai'ly at the larger load. The value of viscos- 
ity used in the calculations reported in table HI was that corresponding to 
the temperature of the lubricant bath. If themial effects become important, 
the isothermal assumption used in deriving equations (15) and (20) is vio- 
lated. But, although the value of the viscosity used in these expressions 
is necessarily somewhat arbitrary, there is evidence from previous studies 
of line or rectangular elastohydrodynamic conjunctions that the film thick- 
ness is determined by the effective viscosity in the inlet region. If this 
viscosity is known with reasonable accuracy, the predicted film thicknesses 
are quite reliable. In addition, at the more severe conditions imposed by 
the larger load the lubricant may no longer behave as a Newtonian liquid. 
This would violate the assumptions used in deriving equations (15) and (20). 

The results of Kunz and Winer (1977) presented in table HI and fig- 
ures 5 and 6 suggest that at large loads there is a more rapid change in 
lubricant film thickness than would be predicted by isothermal elastohydro- 
dynamic theory for elliptical conjunctions as presented in equations (15) 
and (20). This view is supported by the experimental results based on the 
X-ray technique reported by Parker and Kannel (1971) and the results of 
optical interferometry presented by Lee, et al. (1973). This observation, 
however, contradicts the results of Johnson and Roberts (1974), who used the 
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spring dyndtitometer method described in the last section to estimate elasto- 
hydrodynamic film thickness. Johnson and Roberts found that, in spite of 
the approximate nature of their method, the results indicated strongly that 
elastohydrodynamic film thickness as predicted by equations (15) and (i:0) is 
maintaineo up to the highest contact pressures that are practical with ball 
bearing steel. Their results support the more precise measurements of 
Gentle and Cameron (1^73) and extend the range of operating conditions to 
higher rolling speeds, the discrepancy between these investigations in re- 
gard to the influence of load on elastohydrodynamic film thickness has yet 
to be explained. 

Another comparison between experimental findings and theoretical pre- 
dictions can be based on the experimental results provided by Dalmaz and 
Godet (1978). They measured film thickness optically in a pure-sliding, 
circular-contact apparatus for different fire-resistant fluids. An example 
of the good correlation between the theoretical predictions based on the 
formulas developed in the previous section and these experimental results is 
shown in figure 7. The agreement between the experimentally determined 
variation of central film thickness with speed for mineral oil and water- 
glycol lubricants of similar viscosity and the theoretical predictions is 
most encouraging. The ellipticity parameter was unity and the applied load 
was 2.6 newtons for these experiments. Figure 7 also shows that for water- 
glycol the film thickness generated was barely one-third of that developed 
for mineral oil of similar viscosity. This drastic reduction in film thick- 
ness is attributed to the pressure-viscosity coefficient of water-glycol, 
wh-ich is only about one-fifth that of mineral oil. 

Various theoretical predictions are compared with the experimental re- 
sults of Archard and Kirk (1964) in figure 8. These results are presented 
in a form that indicates the influence of side leakage on film thickness in 
elliptical contacts by plotting a film reduction factor r 

against the effective radius ratio Rx/Ry The film reduction factor is 
defined as the ratio of the minimum fi Im -thickness achieved in an elliptical 
contact to that achieved in a rectangular contact formed between two cylin- 
ders in nominal line contact. The theoretical results of Archard and 
CowKing (1965-66), Cheng (1970;, and Hamrock and Dowson (1977), along with 
the theoretical solution for rigid solids by Kapitza (1955), are presented 
in this figure. Archard and Cowking (1965-66) adopted an approach for ellip- 
tical contacts similar to that used by Grubin (1949) for a rectangular- 
contact condition. The Hertzian contact zone is assumed to form a parallel 
film region, and the generation of high pressure in the approach to the 
Hertzian zone is considered. 

Cheng (1970) solved the coupled system of equations separately by first 
calculating the deformations from the Hertz equation. He then applied the 
Reynolds equation to this geometry. He did not consider a change in the 
lubricant density and assumed an exponential law for the viscosity change 
due to pressure. Hamrock and Dowson (1977) developed a procedure for the 
numerical solution of the complete isothermal elastohydrodynamic lubrication 
of elliptical contacts. 

For 0.1 < Rx/Ry < 1 the Archard and Cowking (1965-66) and Hamrock 
and Dowson (1977) predictions are in close agreement, and both overestimate 
the reduction in film thickness evident in Che experimental results of 
Archard and Kirk (1964). For 1 < Rx/Ry < 10 the Hamrock and Dowson 
(1977) predictions of film thickness exceed and gradually diverge from those 
of Archard and Cowking (1965-66). They are also in better agreement with 
the experimental results of Archard and Kirk (1964). Cheng's (1970) theo- 
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retical predictions are optimistic and differ significantly from the experi- 
mental data. Kapitza's (1955) rigid-contact theory is also optimistic about 
the reduction in film thickness compared with the experimental results, but 
it is in closer agreement with experiment than the Cheng (1970) results. 
These comparisons between theory and experiment offer reasonable grounds for 
confidence in the predictions of the complete three-dimensional elastohydro- 
dynamic lubrication theory for elliptical contacts presented in the previous 
section. 


SOFT £HL RESULTS 

The earlier studies of elastohydrodynamic lubrication of conjunctions 
of elliptical form are applied to the particular and interesting situation 
exhibited by materials of low elastic modulus (soft EHL). The procedure 
used in obtaining the soft EHL results is given in Hamrock and Oowson 
(1978). The ellipticity parameter was varied from 1 (a ball-on-plate con- 
figuration) to 12 (a configuration approaching a nominal line or rectangular 
contact). The dimensionless speed and load parameters were varied by one 
order of magnitude. Seventeen different cases used to generate the minimum- 
film-thickness formula are given in table IV. In the table H^-|p corre- 
sponds to the minimum film thickness obtained from applying the EHL ellip- 
tical contact theory given in Hamrock and Oowson (1976) to the soft EHL con- 
tacts. The minimum-f i Irri-thickness formula obtained from a least-squares fit 
of the data was first given in Hamrock and Oowson (1978) and is given here as 


^min * - 0.85 e"^’^^'^) 


In table IV H^^p is the minimum film thickness obtajned from equation 
(21). The percentage difference between H„,^p and Hmin is expressed 
by Oj, given in equation (16). The values of 0^ in table IV are 
within the range -8 to 3 percent. 

It is interesting to compare the equation for materials of low elastic 
modulus (soft EHL, eq. (21)) with the corresponding equation for materials 
of high elastic modulus (hard EHL) given in equation (15). The powers of 
U in equations (21) and (15) are similar, but the power of W is much more 
significant for low-elastic-modulus materials. The expression showing the 
effect of the ellipticity parameter is of exponential form in both equa- 
tions, but with different constants. 

A major difference between equations (21) and (15) is the absence of a 
materials parameter G in the expression for low-elostic-modulus mate- 
rials. There are two reasons for this. One is the negligible effect of 
pressure on the viscosity of the lubricating fluid, and the other is the way 
in which the role of elasticity is simply and automaticol ly incorporated 
into the prediction of conjunction behavior through an increase in the size 
of the Hertzian contact zone corresponding to changes in load. As a check 
on the validity of this, case 9 of table IV was repeated with the material 
changed from nitrile to silicone rubber. The results of this change are 
recorded as case 17 in table IV. 

The dimensionless minimum film thickness calculated from the full 
numerical solution to the elastohydrodynamic contact theory was 
181.8x10“^. and the dimensionless minimum film thickness predicted from 
equation (21) turned out to be 182.5x10-°. This clearly indicates a lack 
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of dependence of the minimum film thickness for low-elastic-modulus mate- 
rials on the materials parameter. 

There is interest in knowing the central film thickness, in addition to 
the minimum film thickness, in elastohydrodynamic contacts. The procedure 
used to obtain an expression for the central film thickness was the same as 
that used to obtain the minimum-film-thickness expression. The central- 
f i litv-thickness formula for low-elastic-modulus materials (soft EHL) as ob- 
tained from Hamrock and Dowson (ly7ti) is 


S * 7.32(1 - 0.72 e-0.i:akj ^,0. 64^-0. 


( 22 ) 


Comparing the central- and minimum-film-thickness equations, (22) and 
(21), reveals only slight differences. The ratio of minimum to central film 
thickness evident in the computed values given in Hamrock and Dowson (1978) 
ranged from 70 to 83 percent, the average being 77 percent. 

Figure 9 shows the variation of the ratio r where 

^*min r is the minimum film thickness for rectangular contacts, with the 
ellipticity pararrieter k for both high- and low-elastic-modulus materials. 
If it is assumed that the minimum film thickness obtained from the elasto- 
hydrodynamic analysis of elliptical contacts can only be obtained to an 
accuracy of 3 percent, we find that the ratio Hmin/Hmin.r approaches 
the limiting value at k « 5 for hich-elastic-modulus materials. For low- 
elastic-modulus materials the ratio approaches the limiting value more slow- 
ly, but it is reasonable to state that the rectangular-contact solution will 
give a very good prediction of the minimum film thickness for conjunctions 
in which k exceeds about il. 

The variation of pressure and film thickness in the direction of roll- 
ing along a line near the midplane of the conjunction is shown in figure 10 

for two values of U. for all the solutions obtained at various speeds, the 

values of the dimensionless load, materials, and ellipticity parameters were 
held fixed. Figure 10 shows that the pressure at any point in the inlet 
region increases as the speed increases. The dominant effect of the dimen- 
sionless speed parameter on the minimum film thickness in elastohydrodynamic 
' contacts for low-elastic-modulus materials (soft EHL) evident in equation 

(21) is reflected in figure 10. Similar results were found for high-elastic- 
modulus materials (hard EHL). 

The variation of pressure and film thickness in the rolling direction 

along a line close to the midplane of the conjunction is shown In figure 11 

for two values of the dimensionless load parameter. Once again the essen- 
tial features of the minimum-f i Im-thickness equation are reflected in this 
figure since a change in the dimensionless load parameter of one order of 
magnitude produces a considerable change in the pressure profile but not 
such a significant change in the film thickness. The small effect of load 
on minimum film thickness is a reasonable and important feature of elasto- 
hydrodynamic lubrication. 

The pressure spikes found when jealing with materials of high elastic 
modulus (hard EHL) are not evident in these solutions for low-elastic- 
modulus materials (soft EHL). The absence of a pressure spike for low- 
elastic-modulus materials is due to the pressures generated for a given load 
in a contact between low-elastic-modulus materials being considerably lower 
than those generated in a contact between high-elastic-modulus materials. 
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CONCLUDING REMmRKS 


The 'lasis of the present lecture was on fully flooded, elastohydro- 
dynamic lubricated elliptical contacts. A fully flooded conjunction Is one 
in which the filin thickness is not significantly changed when the amount of 
lubricant is increased. A brief description of the relevant equations used 
in the elastohydrodynamic lubrication of elliptical contacts was given. The 
most important practical aspect of the elastohydrodynamic theory is the 
determination of the minimum film thickness within the contact. The mainte- 
nance of a fluid film of adequate magnitude is an essential feature of the 
correct operation of lubricated machine elements. The results presented 
show the influence of contact geometry on minimum film thickness as ex- 
pressed by the ellipticity parameter and the dimensionless speed, load, and 
materials parameters. Film thickness equations are developed for materials 
of hi gh elastic modulus, such as metal, and for materials of low elastic 
modulus. Such as rubber. The solutions for high-elastic-modulus materials 
are sometimes referred to as "hard EHL," and the solutions for low-elastic- 
modulus materials as "soft EHL." In addition to the film thickness 
equations that were developed, plots of pressure and film thickness were 
presented. These theoretical solutions for film thickness have all the 
essential features of previously reported experimental observations based on 
optical interferometry. Correlation between theory and experiments was also 
made. 
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APPENDIX - SYMBOLS 


a 


b 


I 


Ol 


£ 


F 

F* 

G 

H 

H 


h 

ho 


P 

P 

R 


r 

S 

u 


u 

U 

Wp 

x,y,X,Y 

a 

n 

'>0 


V 

p 

PO 


semimajor axis of contact ellipse, m 
semiminor axis of contact ellipse, m 

C(Hfnin ~ h^in)/hmin] 100 
modulus of elasticity, N/m^ 


effective elastic modulus 
normal applied load, N 





normal applied load per unit length, N/m 
dimensionless materials parameter, o£' 
dimensionless film thickness, h/RJ^ 

dimensionless film thickness obtained from least-squares fit of 
numerical data 
film thickness, m 
constant, m 
ellipticity parameter 
unit length, m 

dimensionless pressure, p/£' 
pressure, N/m2 
effective radius, m 


radius of curvature, m 

geometric separation, m 

dimensionless speed parameter, nQo/E'Rx 

mean surface velocity in direction of motion, (u^ + ub)/i:, m/s 
dimensionless load parameter, F/£'R^ 

dimensionless load parameter for rectangular contact, F*/£'R- 
coordinate system 

pressure -viscosity coefficient of lubricant, m^^/N 
coefficient of absolute or dynamic viscosity, N s/m^ 

coefficient of absolute or dynamic viscosity at atmospheric 
pressure, N 
Poisson's ratio 
lubricant density, N s^/m^ 
density at atmospheric pressure, N s^/m^ 


Subscripts: 


a 

solid a 

b 

solid b 

c 

central 

min 

minimum 

r 

for rectangular contact 

x.y 

coordinate system 
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TABLE 1. - CrftONOLOGY OF RESEARCH ON ELASTOHYOROOYNAMIC LUBRICATION (EHL) 


OF ELLIPTICAL CONTACTS 


Year 

Major research effort 

Some results 

1974 

(1) Review literature 

(2) Investigate various elasticity 

models 

The elasticity model that assumes a 
uniform pressure over a rectangular 
area was chosen 

1975 

Develop EHL theory with the 
following assumptions: 

(1) Isothermal fluid 

(2) Newtonian fluid 

for the first time a complete theo- 
retical solution was obtained for 
elliptical contacts lubricated 
elastohydrodynamical ly. This in- 
volved successful coupling of the 
elasticity and Reynolds equations 

1976 

Apply EHL theory to hard contacts 
that have 

(1) Fully flooded or 

(2) Starved conjunctions 

Contours describing the pressure 
and film thickness within the 
conjunction were developed 
Film thickness equations were 
formulated for the fully flooded 
and starved situations 

1977 

Apply EHL theory to soft contacts 
that have 

(1) Fully flooded or 

(2) Starved conjunctions 

Same as above 

1978 

Investigate the four fluid-film 
lubrication regimes encountered 
in elliptical contacts: 

(1) Isoviscous rigid 

(2) Viscous rigid 

(3) Viscous elastic (hard EHL) 

(4) Isoviscous elastic (soft 

EHL) 

(1) Film thickness equations were 

formulated for each regime in 
terms of dimensionless geome- 
try, viscosity, and elastic- 
ity parameters 

(2) Maps of the four lubrication 

regimes were obtained 
The results enable both the 
mode of fluid-film lubrication 
and the magnitude of the film 
thickness to be ascertained for 
specified operating conditions 

1979 

& 

1980 

Incorporate the results from the 
previous 5 years of research into 
a book 

EHL results are applied to the 
following applications: 

(1) Ball bearings 

(2) Wheel on a rail 

(3) Gears 

(4) Variable speed drive 

(5) Seals 

(6) Human joints 

(7) Roller bearings 
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T»o. c 1 . _ nuTi FFFECT OF FUlPTlCm. L0«0, SPttO. «N0 MATERIALS PARAMETERS ON MINIMUM FILM THICKNESS FOR 

hard EHL contacts 



- - — " T 

Difference between 
^niin 8(i,in» 

Dl. 

percent 

ReLu t ts 

Case 

1 

El 1 (pticUy 
parameter, 
k 

Dimensionless 
load param- 
eter, 

W 

Dimensionless 
speed -parant- 
eter, 

U 

Dimensi 

mater 

paraa 

G 

on less 

Minimum film thickness | 

lals 

leter. 

Obtained from 
EHL elliptical 
contact theory, 
^min 

Obtained from 
least-squares 

fit. 

Hmin 

1 

d 

3 

4 
b 
6 
7 

U 

y 

10 

1 

l.ifb 

1.5 

1.75 

d 

d.b 

3 

4 
b 
a 

0.110b 

xlO-6 

0.1683) 

10-11 

452 

2 

3.367x10-^ 

4.105 

4.565 

4.y07 

5.255 

5.755 

6.091 

6.636 

6.909 

7.046 

3.514x10-^ 

4.07b 

4.554 

9.955 

5.294 

5.621 

6.196 

6.662 

7.U01 

7.091 

♦4.37 

-.66 

-.24 

♦.96 

♦.74 

♦I.l5 

-1.72 

-.24 

♦.46 

♦.6l 

El 1 ipticity 

11 

id 

1 14 

1 15 
lb 
: 17 

6 

1 


1 

0.2211 

.3beb 

.552a 

.73/1 

.y2l4 

1.10b 

1 . 2 yo 

0.1683x10-11 

1 

45: 

22 

6.492 

6.517 

6.266 

6.156 

6.065 

5.611 

5.657 

6.656 

6.412 

6.225 

6.095 

5.997 

5.916 

5.651 

♦2.53 

♦1.50 

-.69 

-.99 

-1.45 

♦1.84 

♦3.4j 

Load plus 
case 9 

h— 

! lU 

1 

I do 

i 

1 dd 

t 

! k4 

i 

da 

d‘i 

30 

31 

6 

1 

i 

1 

i 


0.7371 


0.08416 

.2525 

.3367 

.420a 

.5ay2 

.6416 

1.263 

1.663 

2.104 

2.525 

2.y4b 

3.367 

4.206 

5.050 

45 

22 

3.926 
6.372 
9.995 
11. bl 
14.39 
16.34 

24.47 
29.75 

34.56 
39.73 

43.47 

47.32 

54.57 

61.32 

3.605 
6.032 
9.769 
11.37 
14.29 
.'6.21 
24.00 
29.16 
33.96 
38.44 
42.69 
4b. 7b 
54.41 
61.59 

-3.06 
-4.0b 
-2.26 
-2.07 
-.69 
-.71 
-1.92 
-1.92 
-1.79 
-3.25 
-1.79 
-1.16 
-.29 
♦ .44 

Speed plus 
case 14 

id 

j3 

34 

— 

6 

6 

6 

0.7216 

.7216 

.24d6 

0.3296 

,y422 

.1122 

2310 

34yl 

6765 

6.931 

17.19 

6.060 

6.936 
17.59 
6. 116 

+0.10 

♦2.33 

1. 

Materials 
plus 
case 9 


TABLE 111. - THEORETICAL AND EXPERIMENTAL FILM THICKNESS 
[From Kunz and Winer (1977). j 


Central film thickness 

Theoretical, 
eq. (20), 

He 

Experi- 
menta 1 , 
He 


speed 

irameter. Theoretical, Experi- Theoretical 
U eq. (20), mental, eq. (15), 

*^min 


Dimensionless load parameter, W, 0.125dxl0~6 



0.1963 x 10-11 

6.84 x 10-^ 

5.7 x 10*8 

" ' ■ '■ 1 

3.87 x 10-8 

2.8 x 10-8 

.3926 

10.9 

9.9 

6.20 

5.7 

.7866 

17.3 

16.0 

9.3 

11.0 

.2637 

11.0 

15.0 

6.58 

8.4 

.5274 

18.9 

22.0 

10.5 

14.0 

1.057 

30.2 

34.0 

17.0 

24.0 

.2268 

8.04 

8.2 

4.53 

5.0 

.4536 

12.8 

12.0 

7.27 

7.5 

.9089 

20.6 

18.0 

11.6 

13.0 

Dimensionless load parameter, U, 

0.9287 x 10-8 


0.1963 x 10-11 

5.96x10“^ 

4.3 x 10-8 

3.33 x 10-8 

2.8 x 10-8 

.3926 

9.50 

7.1 

5.35 

4.3 

.7866 

15.1 

12.0 

8.58 

5.7 

.2637 

10.4 

8.4 

5.68 

2.8 

.5274 

16.6 

12.0 

9.11 

5.6 

1.057 

26.4 

17.0 

14.6 

8.4 

.2258 

7.01 

6.3 

3.92 

2.6 

.4536 

11.2 

9.4 

6.27 

3.8 

,9089 

17.8 

14.0 

10.1 

5.0 


^Lubricant A is poly^lkyl aromatic (o ■ 1.58x10"^ (N/m^)'^; 
no » 0.0255 N s/m^; G ■ 4507); .lubricant B is synthetip 
hydrocarbon (o - 3.UxlO“0 (N/m^)-l; no • 0.0343 N s/m<^; 
G ■ 8874); lubricant C is modified polyphenyl ether 
(o - 1.79x10-8 (N/m^)-l; nfl - 0.0295 N s/m^; G « 5107), 
where a is the pressure-viscosity constant, np is 
the atmospheric viscosity, and G is the dimensionless 
materials parameter, aE'. 
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TABLE iV. - OaTa showing EFFECT OF ELLlPTlClTY, LOaD. SPEED. AND MATERIALS PARAMETERS ON MINIMUM FILM THIUNESS FOR SOFT EHL CONTACTS 


Cose I Ellipiicity I Diniensionless I OinMinslonless Dmiensionless 


Minimum film thickness 


poranieter, 

k 


load paraiit- 
eter, 

w 


speed paranK 
eter, 

U 


0.440bxlO-3 O.iUEBxlO" 


materials 

paraflieter, 

G 


0.4276 


B 

y 

10 

11 

6 

i 

0.4405x10-3 0.05139 

1 .1542 

.2570 

f .05139 

12 

6 

.2202 0.1028 

13 

1 

.6b07 

14 


1.101 

15 


1.542 

16 


2.202 I 

17 

6 

0.1762 0.06169 


Obtained from 
EHL elliptical 
contact theory, 
*^m1n 

BB. 61x10-^ 
i42.5 
170.4 
lBb.7 
20b. 2 
2111.7 
235.2 


0.4276 


0.4276 


Obtained from 
least-squares 
i It. 

Hmin 

yl.OBxlO-^ 
131,2 
160. B 
1B2.4 
20a.B 
224.6 
236.0 


Difference between 
*^m1n “hd Hn,-jn« 

Dl. 

percent 


Results 


Elllpticity 


Speed plus 

case 5 


Load plus 
case 5 


Material 











Dimensionless 



(bl Dimensionless film thickness. 

Figure 4. - Variation ol dimensionless pressure and film 
thickness on X-axis lor three values ot dimensionless 
load parameter. The value of Y is held fixed near axial 
center ol contact. 




Theoretical film thickness, H Theoretical film thickness. 


«xlO'* 



Figure 5, - Theoretical and experimental central and minimum film 
thicknesses for pure sliding. Dimensionless load parameter. W 
aiesaxlO'®. (from Kunz and Wlner, 1977.) 


30x10'* 



Figure 6. ■ Theoretical and experimental central and minimum film 
thicknesses for pure sliding. Dimensionless load parameter, W, 
0. 928x10'®. (From Kunz and Winer, 1977. 1 


Central lilm thickness, 



Figurr 7. - EKect of speed on central lilm thickness 
at constant load (2. 6 Nl for mineral oil and water- 
glycol lubricants of similar viscosity. IFrom Dalmaz 
and Godet, 197&. ) 



Figure 1 - Side- leakage factor for elliptical contacts. 


i 



Ellipttcily p*r«m*t*r. k 


Figure 9. ■ Etiect o» elliplicity parameter on 
ratio dlmensionlesi minimum film thick- 
ness to dimensionless minimum film thick- 
ness lor a line contKt, for EHL high- and 
low-elastic-modulus analyses. 




Figure lOL • Effect of dimensionless speed parameter on dimensionless pressure and film 
thickness on X-a«is. the value of Y is held fixed near axial center of contact. 
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Figure 11. • Edect o( dimensionlesj (o»d perimeter on dimensionless pressure ind lilm 
thickness on X-ixis. Thevilueol Y is held fixed neir ixiilcenter of contict. 




